A non-linear transport method for detecting superconducting stripes 
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We theoretically study the effect of stripe-like superconducting inclusions on the non-linear re- 
sistivity in single crystals. Even when the stripe orientation varies throughout the sample between 
two orthogonal directions due to twinning, we predict that there should be a universal scaling rela- 
tionship between the nonlinear resistivity curves measured at different angles relative to the crystal 
axes. This prediction can be used to verify or rule out the existence of superconducting stripes at 
and above the superconducting transition temperature in cuprate superconductors. 



Study of the high-tcmpcraturc superconductivity 
(HTS) problem has been complicated by proximity of 
the superconducting state to non-superconducting or- 
dered phases. Among them, only antiferromagnetic in- 
sulating state has been unambiguously identified as al- 
ways present in the undoped parent compounds across 
all families of cuprates. The nature of the pseudogap 
regime [l| that emerges upon doping, has remained elu- 
sive, even though there are many theoretical proposals 
for its origin [2|46|. Among the simplest, from the stand 
point of detection, is the possibility that the pseudogap 
regime contains a region of spin and/or charge density 
waves (SDW/CDWV or their strongly-coupled cousin, 
the stripe state 0, |ij . Indeed, neutron scattering com- 
monly detects elastic and inelastic responses characteris- 
tic of static or slowly fluctuating incommensurate SDW 
in majority of cuprate families [3,@1- Charge modulation 
is also observed [9r|l 
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albeit less commonly [12j. Natu- 
rally, charge and spin modulation can locally modify the 
conditions for onset of superconductivity. Moreover, it 
is possible that in the anharmonic regime that charac- 
terizes the stripe state, local pairing mechanisms, which 
are not present in the more common CDW and SDW 
setting, may become operational 13- T^] , in which case 
stripes would be an important contributor to the high 
values of superconducting transition temperature T c . In- 
homogeneous superconductivity could naturally account 
for suppressed superfluid density [3] and for spectral 
weight transfer observed in the optical conductivity of 
underdoped cuprates (Tt| . 

Identification of the stripe states, in contrast to har- 
monic spin or charge modulations has been notoriously 
difficult. Direct experiments that would measure the 
transport anisotropy induced by the rotational symmetry 
breaking caused by stripes are often hindered by sample 
twinning or stripe ordering patterns that restore rota- 
tional symmetry on the measurement length scale [l8j ]. 
In addition, it is believed that in order to be compatible 
with superconductivity, the stripes have to be dynamical, 
i.e., fluctuating on some timescale. 

In this letter we propose a method capable of detecting 
signatures of superconducting stripes diluted in a nor- 
mal matrix despite all these complications. It is based 



on probing the spatial anisotropy of the non-linear trans- 
port, which should be induced by superconducting stripes 
in single crystal samples, even when linear transport is 
completely isotropic. The method takes advantage of the 
fact that unlike the linear response conductivity tensor 
that has to be rotationally invariant in tetragonal sys- 
tems, the non-linear response in general is not. Super- 
conductors have inherently non-linear I-V characteris- 
tics, which makes them ideally suited for the proposed 
method. Therefore, non-linear transport measurements 
can help both to establish the temperature range within 
which superconducting inclusions persist above above T c , 
as well as to determine the local spatial structure of these 
inclusions. 

We qualitatively demonstrate how anisotropic non- 
linear resistivity can arise in a tetragonal system by con- 
sidering a schematic model in Fig. [T^,. Suppose that 
within the Cu02 planes there are domains of supercon- 
ducting stripes, which are primarily oriented along [100] 
and [010] directions. Each superconducting segment has 
a non-linear I-V characteristic, as shown schematically in 
Fig. [TJ). If the characteristic domain size is much smaller 
than the crystal size, the tetragonal symmetry of the sys- 
tem is not broken - the [100] direction is still equivalent 
to [010]. Then, we are immediately forced to conclude 
that linear resistivity must remain isotropic. In partic- 
ular /f? 110 i = P^oo] (superscript "0" indicates that the 
resistivity is taken at zero current). 
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FIG. 1: (a) System with stripes domains aligned along direc- 
tions [100] and [010]. (b) I-V relations used for the super- 
conducting stripes for different values of the critical current 
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Now, let us consider the non-linear resistivity of the 
system. Suppose we apply current j x and measure elec- 
tric field E x , with direction x either parallel to [100] or 
[110] (two high-symmetry directions). For small applied 
current, the superconducting segments have negligible re- 
sistivity, and hence average resistivity of the system the 
smallest. As j x increases, some of the segments quench 
(when the local critical current I c is reached), and the 
overall resistivity increases. We therefore anticipate that 
in both cases, the resistivity will increase non-linearly 
with j x . However, the nonlinear resistivity in two cases 
will not be the same, as can be easily seen from the follow- 
ing argument, which is particularly simple in the dilute 
stripe concentration limit. Suppose that in addition to 
the current in the direction x, there is also a current in 
the y-direction. Qualitatively, it is clear that it should 
make little difference for E x in case x = [100], but a big 
difference when x = [110]. That is because, in the lat- 
ter case, j y current will directly contribute to reaching 
the critical current I c on the SC links involved in the 
transport along x, while in the former, j y primarily sat- 
urates the horizontal superconducting links, which are 
largely irrelevant for the transport along x. Thus we see 
that the non-linear resistivity has to be anisotropic, and 
that the I-V dependence in the coordinate system where 
x = [100] and y = [010] with good accuracy is given by 



E x = p(j x )jx, 
E y = P(jy)jy 



(1) 



That is, the electric field along the stripe pinning direc- 
tion depends only on the current in the same direction. 

In order to obtain the relation between the E and j 
in any other planar coordinate system we only need to 
perform a rotation of Eqs. ([1]) around £-axis by angle 9 
with the matrix TZi, 



E' = ftgp(J)J = K$p(R,i g J')'R.i 6 3' 



where 



p(3) = 



p{J x ) 

P(jy) 



The resistivity tensor in the rotated basis is 

p'(3')=nlp(ni e 3')ni e . 
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Naturally, in the limit of vanishing current we see that 
the resistivity tensors coincide, p' = p, as expected from 
linear response of a tetragonal system. However, at a 
finite current, the relationship is more complex. For in- 
stance, if a current j' is applied along direction x' , the 
electric field response becomes 



E*> = [P(f cos °) cos2 + P(f sin °) sin2 °] f 
E y' = [P(f cos °) ~ P(f sin °)} sin cos 6 j' 



(5) 



In particular, for the high symmetry direction 8 = ir/4 
the resistivity is obtained by dilation of the current axis 
by a y/2 factor, p'(j') = p(j'/\/2). Another notable fea- 
ture is the appearance of the transverse electric field if 
current is applied away from the high-symmetry direc- 
tions, i.e. when 9 is not an integer multiple of 7r/4. The 
induction of transverse electric field in tetragonal systems 
is only possible in the presence of nonlinear response. 

To test the above reasoning, we consider an explicit 
resistor network model that incorporates the described 
qualitative features. The resistors, connecting the near- 
est neighbor sites of a square lattice, iJy, are chosen at 
random to be either "normal" , with constant resistivity 
i?„, or "superconducting", with the current-dependent 
resistivity R S (J S ). The probability of superconducting 
bonds is p. The Kirchhoff's equations 



E 



Vi - Vi 



= 0, 



(6) 



combined with current or voltage boundary conditions 
give a system of nonlinear equations that can be used to 
determine the local voltages {Vi}. In our model calcula- 
tion, for the superconducting resistors we take the I-V 
dependence as 



V S = R S {J S )J S = [R° s 



Rn — R s 
;-4(J,-J ) + I 



J, 



(7) 



which at small currents has resistivity approximately 
< Rn and at large currents J s » I c saturates to 
the resistance of the normal links R n as illustrated in 
Fig. [lb. We will take R n = 1 and R° s = 0.05 unless 
stated otherwise. In the simulations we used a system of 
size 80 x 80 sites with fixed voltage boundary conditions 
in the direction of applied current, and periodic bound- 
ary conditions in the transverse direction. Simulations 
were performed for current applied along the bond and 
diagonal directions. 

In Fig. [2^, we present the result of lattice simulation 
(squares) when current is applied along the bond direc- 
tion for various values of superconducting link concentra- 
tions p and for I c = 2. As expected, for large values of 
applied current J m , the resistivity saturates to the "nor- 
mal" value 1, while at low currents it drops by an amount 
that increases with p. The width of the intermediate non- 
linear region decreases with increasing p. 

Our lattice simulation results can be very accurately 
reproduced by the effective medium theory (EMT). EMT 
is typically applied to study linear random resistor net- 
works, where it well captures, e.g., percolation phenom- 
ena [l9(. Within EMT, one solves exactly the problem 
of a particular random resistor Rq embedded in a per- 
fect matrix composed of identical effective medium resis- 
tors R m . The value of R m is determined sclf-consistently 
from the condition that the voltage drop Vq on resistor 
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FIG. 2: (a) Comparison between EMT and the numerical 
lattice simulation. The lines correspond to Eq. ((9J and the 
squares to the numerical simulation, (b) Verification of the 
"s/2 scaling". Circles and squares correspond to numerical 
lattice simulations for the [100] and [110] directions, respec- 
tively. The dashed lines correspond to the numerical data for 
direction [110] with the current axis divided by y/2. 



i?o, averaged over the distribution, be equal to the av- 
erage voltage drop V m per link of the network. Here we 
apply this procedure to the binary network that contains 
nonlinear resistors. With the average current parallel to 
the bond direction, it leads to the following equations 



Js 

9m{Js) 



2g s (J s ) 



9m{Js) + gs(Js) 



(8) 



(p- \) 2 {9s{J s ) - g n ) 2 +g n g s (Js) 



1/2 



(9) 



where J m is the externally applied uniform current per 
unit cell, J s is the current through a superconducting 
link, which needs to be determined self-consistcntly, and 
for convenience we introduced conductances, gi = I/R4. 
These coupled nonlinear equations can be solved numer- 
ically for any specified form of R S (J). As can be seen 
from Fig. [2^,, non-linear EMT very well approximates 
our lattice simulation results, with the largest deviation 
occurring at the percolation threshold, p = 0.5. 

In Fig. &p we present the comparison of lattice simula- 
tions for the current applied along the bond (circles) and 
diagonal (squares) directions. After rescaling the current 
axis by a factor a/2, the diagonal resistivity (dashed lines) 
matches the resistivity in the bond direction, as was an- 
ticipated from Eq. ([3]). Even at the percolation thresh- 
old, p = 0.5, the scaling works remarkably well. The 
deviation occurs due to the cross-talk between the hor- 
izontal and vertical superconducting bonds, which was 
neglected in Eq. ([1]). We have also performed simulations 
with "longer" superconducting links with lengths of 2 to 
4 lattice constants, which also confirmed a/2 scaling. 

When current is applied in an arbitrary direction rel- 
ative to the crystal axes, in addition to the longitudinal 
resistivity, there is also a finite transverse resistivity. The 
full angular dependence of the longitudinal and trans- 
verse resistivity obtained using Eq. ([5]) for I c = 2 and 



for p = 0.2, 0.8 is presented in Fig. [3l As expected, the 
biggest anisotropy effects are concentrated in the range 
of currents where the resistivity of the superconducting 
links is nonlinear. This transverse response would be 
completely absent in a linear-resistivity of a tetragonal 
system. Therefore, it represents a "null-point measure- 
ment" that can be a very sensitive diagnostic of the pres- 
ence of local rotational symmetry breaking induced by 
superconducting stripes. 
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FIG. 3: Longitudinal (a,c) and transverse (b,d) resistivity as 
a function of applied current (J x , J y ) for p — 0.2 (a,b) and 
p = 0.8 (c,d). 

In order to relate the above model to superconduct- 
ing cuprates we need to know how the parameters 
(p, R®,R n ,I c ) depend on the carrier concentration x and 
temperature T. Qualitatively, we expect that with the 
temperature increasing above T c , the parameters I c and p 
decrease, both vanishing at some temperature T*, while 
R~ and R n increase. Experimentally, T* may correspond 
to the temperature at which the measured linear in-plane 
resistivity experiences a change of slope towards smaller 
values at lower temperatures, the T 2 * °f ^ e ^- poj - If this 
identification is correct, then below T 2 *, a non-linear resis- 
tivity behavior of the kind that we propose should exist. 
We also expect that in the undcrdoped regime, at T c , the 
system is near the percolation threshold, p « 0.5, with 
p > 0.5 at T < T c and p < 0.5 at T > T c . We also 
note that the predicted non-linear transport signatures, 
are distinctly different from another potential source of 
nonlincarity - the Joule heating. Namely, the window 
of currents where non-linear effects due to local super- 
conductivity should exist, shrinks as the temperature in- 
creases above T c , vanishing at some T* , while Joule heat- 
ing should increase with increasing current indefinitely 
and should be sensitive to the details of the thermal cou- 
pling to environment. Moreover, Joule heating is not 
expected to cause any spatial anisotropy in resistivity 
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Another issue that we can comment on within our su- 
perconducting percolation model is the apparent discrep- 
ancy between the superconductin g flu ctuation regime, as 
extracted from the bulk probes 1211 and from the AC 
transport measurements [17|, [22|, |23[ . In the AC trans- 
port, one attempts to detect features characteristic of 
superconductivity, such as the superfluid density, and 
tracks their disappearance as a function of increasing 
temperature. Experiments consistently find that super- 
conducting fluctuation range above T c , extracted this 
way, is much more narrow than the one obtained from the 
bulk measurements, such as diamagnetic response plj ]. 
Within the uniform fluctuating superconductor scenario 
these differences arc difficult to reconcile (23j| . However, if 
we assume that the superconductor is intrinsically inho- 
mogeneous, as in the model considered above, the rapid 
disappearance of the transport signatures of supercon- 
ductivity can be interpreted in terms of reduction of the 
superconducting fraction p below the percolation thresh- 
old. Indeed, let us assume that the conductivity of the 
normal links is g n = ar n and of the superconducting 
links is g s ~ ar s /(l + iu>t s ) in the relevant frequency 
range, ujt h <C 1, where r s is the relaxation time in the 
superconducting links and r„ is the normal state relax- 
ation rate (t s 3> t„). The effective medium conductivity 
g m can be obtained from EMT, Eq. ((9]). We find, above 
percolation threshold, g m ~ (2p — l)g s , while below - 
5m ~ — 2p), the latter having only a very small 

imaginary part. The crossover occurs in a narrow range 
of Sp ~ y/ijJTn 1, i.e., the superconducting contribu- 
tion disappears very rapidly as the fraction of the su- 
perconducting links goes below the percolation threshold 
p = 1/2 near T c , even though disconnected superconduct- 
ing inclusions may still persist to much higher tempera- 
tures and contribute to other superconductivity-sensitive 
measures. In Fig. 2] we show schematically the behavior 
of j m (w) expected from the EMT, which agrees with the 
trends observed experimentally. 
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FIG. 4: Frequency dependence of (a) Real and (b) Imaginary 
parts of conductivity obtained from effective medium theory 
for systems with different stripe concentrations p. t„/t s — 
100. 



The lattice percolation model presented here, while 
crude, provides a transport-based method to test the ex- 
istence of local superconducting inclusions - supercon- 
ducting stripes - in the cuprates at and above T c . In 
combination with the bulk probes, which are sensitive to 
the local superfluid density (oc I c ) and the volume frac- 
tion p, it may help to shed light on the nature of the 
pscudogap regime, from which the high-temperature su- 
perconductivity emerges. 
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